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An axial next-nearest-neighbor Ising (ANNNI) model is studied by using the non-equilibrium 
relaxation method. We find that the incommensurate stripe phase between the ordered phase and 
the paramagnetic phase is negligibly narrow or may vanish in the thermodynamic limit. The phase 
transition is the second-order transition if approached from the ordered phase, and it is of the 
Kosterlitz-Thouless type if approached from the paramagnetic phase. Both transition temperatures 
coincide with each other within the numerical errors. The incommensurate phase which has been 
observed previously is a paramagnetic phase with a very long correlation length (typically £ > 500). 
We could resolve this phase by treating very large systems (~ 6400 x 6400), which is first made 
possible by employing the present method. 
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I. INTRODUCTION 

Of late years, incommensurate (IC) stripe structures 
have been interesting subjects in various physical phe- 
nomena. As the typical examples, we may list an al- 
loy ErgoYio- x La x which contains heavy rare earth met- 
als, Er, Tm, [fflQ the incommensurate phase of dielectric 
material such as Pb(Zri_ x Ti x )03 and NaNC>2, f^Q] and 
the stripe structure in Cu02 planes of oxide supercon- 
ductors. H In Erg Yio_ x La x , the longitudinal incom- 
mensurate oscillatory phase appears between the para- 
magnetic phase and the ordered phase. The aligned 
holes (domain walls) separate antifcrromagnetic stripes 
in CuC>2 planes of oxide superconductors, then the spin 
and charge arc modulated. In such systems, coop- 
erative effects of fluctuation and frustration are con- 
sidered to play important roles. Thus, we sometimes 
treat them with the axial next-nearest-neighbor Ising 
(ANNNI) model as the simplified theoretical model. For 
instance, when a uniaxial anisotropy is strong in the di- 
electric material, the Hamiltonian is equivalent to the 
ANNNI model if we only consider the dipole interactions 
up to the next-nearest-neighbor distance. The phase di- 
agram of Pb(Zri_ x Ti x )C>3 obtained by experiments ||] 
agrees with that of the three-dimensional ANNNI model 
obtained by the mean-field approximation. || In cop- 
per oxide materials, CasLa6Cu2404i and Ca2Y2CusOio, 
the Cu-O-Cu chains with ferromagnetic nearest-neighbor 
and antiferromagnetic next-nearest-neighbor interactions 
are aligned on two-dimensional planes whose interchain 
interactions are antiferromagnetic. Furthermore, the 
spins on this plane are predicted to have a strong Ising 
anisotropy. 0,|l) Consequently, we may treat these cop- 
per oxide planes as the two-dimensional ANNNI model, 
between which the conduction electron planes exist. 

In the ANNNI model, there are exchange interac- 
tions up to the next-nearest-neighbor pairs along one 
axis, while there are only the nearest-neighbor inter- 
actions along the other axes. Most commonly, we 
take a convention that the nearest-neighbor interactions 
are ferromagnetic and the next-nearest-neighbor ones 
are antiferromagnetic, which cause frustration. When 
the next-nearest-neighbor interactions (frustration) are 
weak, a ferromagnetic state is the ground state, and a 
paramagnetic-ferromagnetic phase transition occurs. On 
the other hand, when frustration is strong, the ground 
state is the antiphase state (abbreviated by (2)), which is 
a commensurate(C) stripe structure of two up-spins and 
two down- spins like Ttlltt- It is widely accepted that the 
incommensurate stripe phase exists between the param- 
agnetic phase and the antiphase, and that the succes- 
sive phase transitions (paramagnetic-IC-antiphase) take 
place; there is a "finite" incommensurate phase where the 
up-spins and down-spins are aligned with a period longer 
than two. 



Although only a topology of the phase diagram is 
known by the mean-field approximation in the three di- 
mensional model, l| estimates for the phase transition 
temperatures in two dimensions have been done by sev- 
eral approximation theories [p~2|— p~7[ and numerical simu- 
lations. 1 1 8 20| However, values of the phase transition 
temperatures scatter much depending on the methods 
employed. One of the reasons that make the estimate dif- 
ficult is a lack of a reliable numerical simulation. Because 
of frustration, the correlation time of Monte Carlo (MC) 
simulations becomes very long near the critical point in 
large systems. Thus, we are hardly able to reach thermal 
equilibrium states within reasonable time steps. A pow- 
erful method called a cluster heat bath (CHB) method 
ifiof was developed to reduce the correlation time, how- 
ever, the system size accessible by this method is re- 
stricted up to about 64 x 128 lattice sites. 

Recently, a new method using a MC simulation has 
been developed: the non-equilibrium relaxation (NER) 
method. |^Hp2| We are able to understand phase transi- 
tions from the differences of behaviors in non-equilibrium 
relaxation processes, which have been discarded in con- 
ventional MC simulations. Since we do not wait until 
the equilibrium is realized, we can use the CPU time to 
enlarge the system sizes. Therefore, it makes possible 
to treat large systems that can not be possible by other 
methods. Accordingly, there expected to be little finite- 
size effect in the obtained data, and thus we can regard 
them as those of the infinite systems. By reasons stated 
above, the NER method can be effective especially in 
systems with a slow dynamics which a very long equili- 
bration is necessary. p3[-^H 

In this paper, we study the two-dimensional ANNNI 
model by the NER method in order to determine the 
successive phase transition temperatures. However, the 
obtained results suggest that the IC phase may disap- 
pear within limits of numerical errors. Moreover, we 
found a very exotic phase transition, which seems to be 
the Kosterlitz-Thouless (KT) transition |33) if it is ob- 
served from the high temperature side, and seems to be 
the second-order transition if observed from the low tem- 
perature side. This evidence suggests that the frustration 
parameter does not serve as an asymmetric parameter 
p4]~^6| which explicitly favors the IC structure. 

We describe the model Hamiltonian and the NER 
method in Sec. D, and present simulational results in 
Sec. III. Section fvl is devoted to conclusions. 



II. MODEL AND METHOD 

A. Two-dimensional ANNNI model 

The two-dimensional ANNNI model is described by the 
following Hamiltonian, 
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where Jo(> 0) is the nearest-neighbor interaction along 
the chain direction, which is the direction that has no 
frustration. Ji(> 0) and J 2 (< 0) are respectively the 
nearest-neighbor and the next-nearest-neighbor interac- 
tions along a direction perpendicular to the chain direc- 
tion (axial direction that has frustration), and S x . y = ±1. 
In this paper, we fix J = J\ for simplicity, and impose 
the open boundary conditions along the axial (y) direc- 
tion, while we use the periodic boundary conditions along 
the chain (x) direction. We define a ratio between the 
nearest and the next-nearest-neighbor interactions along 
the axial direction as k(= —Ji/Jx). Our interest is re- 
stricted to the region k > 1/2, where the successive phase 
transitions (antiphase — > IC — > paramagnetic) has been 
considered to occur. 

The fermion approximations assumed that the system 
has straight domain walls along the chain direction at the 
intermediate temperatures. [ p"3| That is, domains whose 
periods are longer than two appear among the commen- 
surate antiphase domains. In this case, the spin structure 
becomes incommensurate along the axial direction, while 
the spins along the chain direction order ferromagneti- 
cally. On the other hand, it was postulated that domain 
walls run along the axial direction in the interface free 
energy method of Muller-Hartmann and Zittartz(MHZ). 
[ p2| , pT| Namely, the correlation of spins ordered ferromag- 
netically along the chain direction is destroyed by these 
domain walls. 

Sato and Matsubara |^o| discussed using the CHB sim- 
ulation that the transition temperature (T c2 ) between the 
antiphase and the IC phase agrees with that of the free- 
fcrmion approximation while the transition temperature 
(T c i) between the IC phase and the paramagnetic phase 
is close to T c2 obtained by MHZ. Thus, they considered 
that the domain walls penetrate into the system along 
the chain direction at T c2 as the temperature increases 
from the ground state. The ferromagnetic correlations 
along the chain direction remain until they are destroyed 
at T c i and the paramagnetic state is realized. This idea 
means that the IC phase exists between the penetration 
of incommensurate domains in the axial direction and 
the disappearance of the ferromagnetic chain correlation. 
Therefore, the spin structure along the axial direction is 
incommensurate, while it is ferromagnetic along the chain 
direction in the IC phase. It has been considered that the 
phase transition between the paramagnetic phase and the 
IC phase is the KT transition and that between the IC 
phase and the antiphase is the second-order transition. 

We consider following two quantities to clarify the 
phase transitions. One is the antiphase magnetization 
defined by 



,r,y 



where N is the total number of spins in the system, 
Sx,y{t) denotes the spin value of the site (x,y) at time 
t, and S X:V ((2}) represents the antiphase ordered state. 
The antiphase magnetization takes a finite value in the 
antiphase, but vanishes in the IC phase and the paramag- 
netic phase, since the domain walls destroy the antiphase 
state at the lower transition temperature, T c2 . Therefore, 
this parameter is employed to estimate T c2 . 

Here, we should note that the antiphase magnetiza- 
tion is not a relevant order parameter that decays alge- 
braically at the transition temperature. Since the ele- 
mentary excitation is a domain wall which readily per- 
colates the chain direction, the antiphase magnetization 
decays exponentially once the domain wall penetrate into 
the system. Usually, the density of the domain walls has 
been used as an order parameter. This parameter is con- 
sidered to diverge algebraically at the transition temper- 
ature as (T — T c2 )p , which is equivalent to that the cor- 
relation length diverges algebraically as £ ~ (T — T c2 )~ v . 
In this point, the phase transition between the IC phase 
and the antiphase is the second-order transition. In the 
Monte Carlo simulations, the correlation length is re- 
lated to the characteristic time as r ~ £ 2 , where z is 
called the dynamic exponent. Therefore, the character- 
istic time also diverges at T c2 as r ~ (T — T c2 )~ zu . We 
can extract it from the relaxation of the antiphase mag- 
netization: the time that the antiphase magnetization 
begins to decay exponentially is the characteristic time 
that the domain wall penetrates into the system. We ex- 
ecute the finite-time scaling to obtain the characteristic 
time at each temperature and estimate T c2 as its diverg- 
ing temperature. A concrete procedure of the finite-time 
scaling is explained in Sec. II D. The antiphase magne- 



tization is an extensive variable so that it shows better 
accuracy as the system size is enlarged because of the self 
averaging. Thus, we use the antiphase magnetization to 
estimate the lower transition temperature, T c2 . 

The other quantity is the layer magnetization defined 



as 



mi(t) 



(3) 



where L x and L y are the length of the system along the 
chain and the axial direction, respectively. The layer 
magnetization reflects the spin order along the chain di- 
rection. In the IC phase, the incommensuration is real- 
ized in the axial direction. The spin arrangement along 
the chain direction should order ferromagnetically or at 
least be critical. The free energy of this state, which is of 
the order of — L 1+ ~ , is always lower than that of the state 
with the spins disordered only along the chain direction, 
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which is of the order of —L, if v is positive. Accord- 
ingly, the layer magnetization vanishes exponentially in 
the paramagnetic phase, but takes a finite value or be- 
haves in a power law in the IC phase. Therefore, this 
order parameter is used to estimate the upper transition 
temperature, T c %. We determine the transition tempera- 
tures by using these two quantities in this paper. 



B. Non-equilibrium relaxation (NER) method 

Phase transitions occur in the infinite-size limit(L — > 
oo) and in the equilibrium limit (t — > oo). Because we 
cannot take both limits at the same time in the simula- 
tion, the equilibrium limit has been taken first conven- 
tionally, and then, we take the infinite-size limit by using 
the finite-size scaling. However, the dynamics of the sim- 
ulations are very slow in the frustrated systems, which 
causes a very large correlation time. From a time-space 
relation, r ~ £ z , this means that the correlation length 
is also very large. As will be mentioned in Sec. [II, we 
estimate the correlation length, £ x ~ 500, for k = 0.8 
and T = 1.40, which resides in the IC phase of previ- 
ous phase diagram. In this situation, reliability of the 



finite-size scaling might become doubtful. In this paper, 
we follow a completely alternative approach to the ther- 
modynamic limit, i.e., we observe the relaxation of the 
infinite-size system to the equilibrium state. In order to 
extract the equilibrium properties, the finite-time scaling 
analysis is utilized instead of the conventional finite-size 
scaling. This approach is known as the non-equilibrium 
relaxation (NER) method. 
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31] Actually, we pre- 
pare a very large lattice and observe the relaxation of 
physical quantities. The simulations are stopped before 
the finite-size effect appears. Accordingly we can regard 
the systems as the infinite systems. 

Using the NER method, we can estimate the phase 
transition temperature and the critical exponents by ex- 
amining behaviors of the relaxation processes to the ther- 
mal equilibrium state (non-equilibrium relaxation pro- 
cesses). The analysis is based on the dynamic finite-size 
scaling hypothesis of the free energy that 



F(e, h, L, t) = L- d F(eLi , hL d ~» , tL~ z ), 



(4) 



where e(= (T — T c )/T c ), h, L, t are the relative temper- 
ature, the symmetry breaking field, the system size and 
time, respectively, v and /3 denote static exponents, while 
z is a dynamic exponent, and d is a dimension of space. 
|57j The order parameter is given by a derivative of the 
free energy with the field: 



m(t) 



dF 
~dh 



L--F(eL-,tL- z ) 



= t-^F(st~) 
ex t~~^ (e = 0, L 



where we have set L ~ t x l z at the transition point (e = 0) 
because the characteristic time scale r and the correlation 
length £ should scale as £ ~ t 1 / 2 . Equation (H) describes 
that the order parameter decays in a power law with time 
at the transition temperature. On the other hand, the 
order parameter behaves exponentially at temperatures 
away from the critical point. Thus, we are able to es- 
timate T c by this difference. Actually, we measure the 
order parameter m(t) at each MC step (i) started from 
an ordered spin configuration at a given temperature. 
We repeat this MC run by changing the random number 
seeds, and m(t) is averaged over these independent MC 
runs. The temperature at which the m(t) curve decays 
in a power law is the transition temperature. 

We use the local exponent to ascertain whether the 
order parameter decays in a power law or exponentially, 
which is defined by 



A(t) = 



d log m(i) 



dlogf 



(6) 



When we plot the local exponent against 1/t, it diverges 
to infinity for T > T c , it converges to zero for T <T C and 
to a finite value (^ 0) at T = T c in the limit of 1/t — » 0. 
The upper bound of T c is the lowest temperature that 
A(t) diverges, and the lower bound is the highest tem- 
perature that A(t) decays to zero. The convergent value 
of the local exponent at T = T c is the critical exponent 
A = 0/zv from Eq. (|). 



C. NER of fluctuation 

Here, we describe the NER of fluctuation. f^,|||||,||] 
The susceptibility is written by differentiating Eq. (0) 
with the symmetry breaking field twice, 



d 2 F 



dh 2 



oc (m(t) 2 ) - {m(t)f 



h=0 



cx t~~~ — t~ , 

where we have used the scaling relations: 

dv = 2 — a 
a + 2(3 + j = 2. 



(7) 



(8) 



oo), 



(5) 



The susceptibility diverges in a power law at the tran- 
sition temperature. Thus, we are able to estimate the 
transition temperature and the critical exponents from 
the susceptibility. It is also noticed that the NER of fluc- 
tuation does not require us to start with a symmetry- 
broken ordered state. The quantity of fluctuation always 
takes a definite value and diverges at T c even though the 
symmetry is not spontaneously broken. Therefore, we 
can start from a paramagnetic state or any state in this 
scheme. This is especially useful when an ordered state is 
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not known yet, or it is difficult to realize. The layer mag- 
netization defined by Eq. (^) is equivalent to the second 
derivative with respect to a local field along a single chain 
in the paramagnetic phase. When the simulation starts 
from the antiphase, the first derivative term remains fi- 
nite in the NER process even though the temperature is 
in the paramagnetic phase. Therefore, the layer magne- 
tization is not regarded as the susceptibility. The NER 
of the layer magnetization is within a scheme of the NER 
of fluctuation, only when the simulation is started from 
the paramagnetic state. 



D. Finite-time scaling 

In case when the NER function does not begin to decay 
algebraically within a reasonable time, it is difficult or 
almost impossible to estimate the transition temperature 
directly by the local exponent. Even in such a situation, 
we are able to determine it by using the finite-time scaling 
analysis, which is a direct interpretation of the finite- 
"size" scaling by a relation £ ~ r 1 / 2 . We present the 
finite-time scaling relation |3Q,pl| as follows, 



m(e,t) = t- x m(t/T(e)), 



(9) 



where A = (3/zv, and r(e) denotes a relaxation time at 
the relative temperature e. Since the relaxation time di- 
verges algebraically in the case of the second-order phase 
transitions, the relation between the relaxation time and 
the relative temperature is described by 



r(e) = Ae 



(10) 



On the other hand, in the case of the KT transition, |33j ] 
it is considered that the relaxation time diverges expo- 
nentially, which we assume that 



r(e) = Aexp(B/Vi). 



(11) 



Now, we describe how we actually estimate the tran- 
sition temperature by using Eq. We use only 
data which is clearly in the paramagnetic phase; T > T c . 
First, we plot m(t)t x at various temperature against 
t/r(s) by using Eq. (||) to determine A and r(e) so that 
all data points fall on a single curve. Next, we plot r(e) 
against e, and fit the points to a smooth curve as we 
change T c and zv by using Eq. ( [To|) for the second-order 
transition. In the KT-transition case, we use Eq. ( [TT] ) in- 
stead. The temperature at which the least-square fitting 
error becomes minimum is the most probable estimate 
for T c . The phase transition between the paramagnetic 
phase and the IC phase is the KT transition and that 
between the IC phase and the antiphase is the second- 
order transition in the two-dimensional ANNNI model. 
Therefore, we use Eq. (^) for the upper transition tem- 
perature, and Eq. ( |To| ) for the lower one. 



III. SIMULATION AND RESULTS 
A. NER from the antiphase state 

We examine the phase transition temperature for k — 
0.6 and 0.8. The temperature is measured in a unit of 
Ji . Here, the NER of two quantities are presented by the 
simulation started from the antiphase state. 

First, we observe the NER of the layer magnetization 
to check the relation between the finite size, L x , and 
the correlation length, £ x , along the chain direction. We 
change the size along the chain direction from L x = 799 
to L x = 25599, while the length along the axial direc- 
tion is fixed to L y = 800. Figurel shows the NER of the 
layer magnetization, Eq. (||), for k = 0.8 and T = 1.40. 
This temperature belongs to the IC phase (KT phase) in 
previous investigations as summarized in Table Q. When 
the system size is small (N — 799 x 800), the relaxation 
roughly looks like a power-law decay, which misled us to 
the KT phase. As the system becomes larger, the relax- 
ation exhibits an exponential decay confirming us that 
the system is in the paramagnetic phase. Note that the 
convergence to a finite value is due to the finite-size effect. 
By its definition, the equilibrium value of the layer mag- 
netization is roughly estimated as (, X /L X in the param- 
agnetic phase. As shown in Fig.l, the convergence value 
takes a half-value if the system size is doubled. Thus, 
we can estimate the correlation length £ x ~ 500 in this 
system. From this figure, we can also understand the re- 
lation between the finite size-effect and the time-effect; 
the effective time that the system behaves as the infinite 
size. For example, it is about 15000 MCS in the system 
with L x = 6399. After this time scale, finite size effect 
appears in dynamics. In this subsection, we use the lat- 
tice of N = 6399 x 6400 and determine the observing time 
at each temperature until which the relaxation curve do 
not begin to bend to an equilibrium value. We average 
24 ~ 32 independent MC runs at each temperature. 

Figure2 is a raw data of the NER of the antiphase mag- 
netization, Eq. (0), for k = 0.6. At all the temperatures, 
the antiphase magnetization clearly decays exponentially, 
which guarantees that the transition temperature must 
be lower than 0.98. As mentioned in Sec. II A, the an- 



tiphase magnetization is not a relevant order parameter, 
and so it decays exponentially as soon as domain walls 
penetrate into the system. Therefore, we estimate the 
transition temperature by using the finite-time scaling. 
At first, we determine the exponent, A, and the relaxation 
time, T(2) (e), at each temperature so that the scaled data, 
rri(2)(t)t x , fall on a single scaling function when plotted 
against t/rio\{e) (Fig. 3(a)). Here, the relaxation times 
are normalized by the value of T = 1.06, and are listed 
in Table |l[ The exponent, A, that gives best scaling is 
A = 0.015. Excellence of the scaling shown in Fig. 3 may 
yield the validity of the finite-time scaling hypothesis. 
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Next, we estimate T c2 by Eq. (|T(]) because it is pre- 
dicted that the phase transition between the antiphase 
and the IC phase is the second-order. We plot t^)(s) 
against e(= T ^ c2 ) as changing T c2 and find a T c2 that 
gives the best linearity in the log-log scale as shown in 
Fig.3(b). We obtain T c2 = 0.89 ± 0.02 and zv = 4.07. 
The error ±0.02 is a range of temperature in which the 
data points clearly fall on the fitting line. The exponent 
has a range within which we can perform a good scaling 
as A = 0.000 ~ 0.030. The transition temperature takes 
the same value irrespective of the choice of A. Evidences 
such as the algebraic divergence of r^ 2 ) (e) and excellence 
of the finite-time scaling support that this transition is 
the second-order. This is a clear distinction from the 
three dimensional model, where the lower transition is 
considered as the first-order. 

We estimate the upper transition temperature, T c \, us- 
ing the layer magnetization, Eq. (|^), in the same way as 
mentioned above. Figure4 shows raw data of the NER 
of the layer magnetization, mi(t), for k = 0.6. Here, 
the layer magnetization clearly decays exponentially, and 
thus the T c2 must be less than 0.98. The finite-time 
scaling is shown in Fig. 5(a). The obtained relaxation 
time, r;(e), is also presented in Table. |j[ Since it is pre- 
dicted that the phase transition between the IC phase 
and the paramagnetic phase is the KT transition, we fit- 
ted Ti(e) by Eq. ( [Ti"| ) in Fig. 5(b), by which we obtain 
T c i = 0.89 ± 0.02 and B = 3.44. If we assumed that the 
phase transition is the second-order transition, and per- 
formed the fitting by Eq. ([l0|), the estimated T c \ becomes 
far from a physically meaningful value. In consequence, 
we confirm that the upper phase transition is of the KT 
type. Actually, Sato and Matsubara Q performed the 
finite-size scaling of the layer magnetization supposing 



m,xL' = Y\L- X cxp(S£-°- 5 )], 



(12) 



which gave T c \ = 1.16 and r\ — 0.25 for k = 0.6. 
This scaling form is equivalent to the finite-time scaling, 
Eq. (g) and (0), if we admit L = t 1 ^ and A = jj/z. The 
difference of the obtained transition temperatures can be 
attributed to the difference of the system sizes. We esti- 
mate A = 0.015 by minimizing the normalized residual, 
however, this value does not correspond to (2 — rf)/z. 
Because we start simulations from the antiphase state, 
the layer magnetization contains contributions from the 
antiphase magnetization, (to^ 2 )) 2 - 

We have obtained the T c2 and the T c \ for k — 0.8 in 
the same way. In the finite-time scaling plot, we used the 
data at five temperature ranging T = 1.40 ~ 1.48, and 
obtained T c2 = 1.32 ± 0.02 and T cl = 1.31 ± 0.02. 

The phase transition temperatures, T c \ and T c2l coin- 
cide within limits of error both for k = 0.6 and 0.8. This 
means the IC phase may vanish in the thermodynamic 
limit. Because this is a rather daring conclusion, we must 
confirm it from another point of view. Actually, it might 



be dangerous to estimate the transition temperature be- 
tween the paramagnetic phase and the IC phase started 
from the ground state. 



B. NER from the paramagnetic state 

We start the simulation from the paramagnetic state 
and observe the NER of fluctuation; the layer magneti- 
zation. Since the observable is a quantity of fluctuation, 
it is necessary to take much more sample averages com- 
pared to the NER of antiphase magnetization. There- 
fore, the system size is restricted to N = 1999 x 2000 
at most. The NER function of fluctuation diverges alge- 
braically at T = T c , diverges exponentially at T < T c and 
remains finite at T > T c . We find the phase transition 
temperature from these differences. It is noticed that the 
NER function finally converges to a finite value because 
the system is finite. Therefore, we must be aware of the 
range of time not to observe the finite-size effect. This 
is a time scale that the correlation length reaches the fi- 
nite system size. We compare the NER functions of two 
different sizes (1999 x 2000 and 1599 x 1600/1599 x 1600 
and 999 x 1000/1599 x 1600 and 799 x 800 et al.) for 
every data point, and estimate this crossover time until 
which two curves fall on the same line and the finite-size 
effect does not appear. An example of this comparison is 
shown in Fig. 6. The NER functions of the layer magne- 
tization from the paramagnetic phase, mi(t), are plotted 
for three sizes, N = 1999 x 2000, N = 1599 x 1600 and 
N = 1199 x 1200 at T = 0.92 and K = 0.6. Three curves 
fall on the same line until 25000 MCS. After this crossover 
time, however, the curve of N — 1199 x 1200 bends down 
from the other two curves, probably because the corre- 
lation length reaches 1200 at this temperature. So, we 
have to discard the data of N = 1199 x 1200 after 25000 
MCS. Comparing two curves of N = 1999 x 2000 and 
N = 1599 x 1600, we can use the data of N = 1599 x 1600 
until at least 10 5 MCS. Table III shows the crossover 
MCS at each temperature for k — 0.6 and 0.8. Thus, we 



use a system of a proper size for a proper observing time 
at each temperature. We take averages over 1000 ~ 7000 
independent MC runs. 

The NER of the layer magnetization, m/(i), from the 
paramagnetic state are shown in Fig. 7 for (a)K = 0.6 
and (b)re = 0.8. Figure8 shows the corresponding local 
exponents. In Fig. 8(a), the exponent decreases for T > 
0.92 and diverges for T < 0.88. At T = 0.90, it converges 
to a finite value. In consequence, we predict that the 
phase transition temperature between the paramagnetic 
phase and the IC phase is T cl = 0.90 ±0.02. For k = 0.8, 
we obtained T c \ = 1.325 ± 0.025 as shown in Fig. 8(b). 

Furthermore, we also analyze the upper transition tem- 
perature, T c i, using the finite-time scaling as shown 
in Figs. 9 and 10. We show the scaling in Fig. 9 for 
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(a)« = 0.6 and (b)/c = 0.8. Since the layer magne- 
tization diverges algebraically as i — ~ at the transition 
temperature, we plot mi(t)t~ x against t/Ti{e). All the 
curves excellently fall on the same line, by which we ob- 
tain the relaxation time as shown in Fig. 10. Here, we 
adopt A = (2 — rj)/z = 0.49 that is a value which the local 
exponent converges to (see Fig. 8). If we admit the KT- 
criterion rj = 1/4, the dynamic exponent is estimated as 
z ~ 3.6. Next, we fit the relaxation time using Eq. (11). 
FigurelO shows the best least-square fitting for (a)« = 0.6 
and (b)n = 0.8. We obtain the transition temperature 
T cl = 0.890±0.015 for k = 0.6 and T cl = 1.300±0.013 for 
K = 0.8. These transition temperatures are again consis- 
tent with the values estimated from the local exponents 
of the NER from the paramagnetic state, and those by 
the scaling from the ground state. 

Therefore, we are able to conclude that T c \ and T c2 
coincide with each other within limits of error in both 
cases of k — 0.6 and 0.8. The phase transition tempera- 
ture between the paramagnetic phase and the IC phase 
is T c i = 0.895 ± 0.025 and that between the antiphase 
and the IC phase is T c2 = 0.89 ± 0.02 for k = 0.6, and 
T cl = 1.32 ± 0.03 and T c2 = 1.32 ± 0.02 for k = 0.8. 
Since these temperatures are very close to each other, it 
is suggested that the IC phase does not exist or it is very 
narrow even if it exists. We need to pay attention to a 
fact that Monte Carlo simulation is not able to exclude 
an very tiny temperature region. 

IV. CONCLUSIONS 



It has been considered that the successive phase tran- 
sitions with a finite IC phase take place for k > 0.5 in 
two-dimensional ANNNI model. In this paper, we es- 
timated the phase transition temperatures by applying 
the NER method, and found that the T c \ is equal to T c2 
within limits of errors. This is a very exotic phase transi- 
tion, which is the KT type if approached from the high- 
temperature side, and is the second-order if approached 
from the low-temperature side. Therefore, we specu- 
late successive phase transitions with an infinitesimal IC 
phase may occur in this system. 

In the studies of the C-IC transitions in two- 
dimensional systems, it has been investigated the systems 
with an asymmetric parameter which explicitly favors the 
IC structure. |34|-|3(| For a finite value of the parameter, 
there exists the finite IC phase between the commensu- 
rate phase and the paramagnetic phase. In the limit of 
vanishing the parameter, the IC phase shrinks to the in- 
finitesimal. The frustration parameter, k, of the ANNNI 
model has been considered as the asymmetric parame- 
ter based on the approximate theory ||T^— |l7|l valid only 
at low temperatures, and the small-scaled Monte Carlo 
simulations. lia-GG] However, what is actually favored by 



frustration is the creation of the domain wall. Between 
the creation of the domain wall and the realization of the 
IC structure, there are many conditions to satisfy, which 
have been supposed in the free-fermion approximation. 
p3[ One of these is the spin correlation along the ferro- 
magnetic chain. In the two-dimensional ANNNI model, 
the non-frustrated direction is only one dimension, and 
thus the spin correlation along the ferromagnetic chain 
can be easily destroyed. Therefore, we question regarding 
the frustration parameter as the asymmetric parameter. 
If we consider that these two are not related with each 
other and the asymmetric parameter is zero in the present 
model, the width of the IC phase becomes infinitesimal. 
This is what we observed in this paper. 

Here, we show the comparison of the obtained tran- 
sition temperatures with the previous ones in Table Q. 
It is recognized that our T c \ is lower than any other 
ones, though T c2 is consistent with each other. This 
can be explained by the difference of the finite-size ef- 
fect. The phase transition between the IC phase and the 
paramagnetic phase is confirmed to be the KT transition 
while the one between the antiphase and the IC phase 
is the second-order. The correlation length diverges al- 
gebraically against the temperature in the latter, while 
it diverges exponentially in the former case. If a system 
size is small compared with the correlation length and 
the accuracy of the numerical data is not enough, even 
the finite-size scaling analysis may mislead to a wrong T c , 
where the finite but very large correlation length reaches 
the finite system size. Therefore, it is very likely that the 
KT transition temperature obtained previously is over- 
estimated. We may rather easily obtain the phase tran- 
sition temperature accurately in the second-order tran- 
sition, because the divergence is algebraically. In the 
two-dimensional ANNNI model, the correlation length 
along the chain direction is £ x ~ 500 at T = 1.40 and 
k = 0.8 which is the temperature a little higher than the 
KT transition temperature and the correlation length is 
far larger than the system sizes ever simulated previously 
(N = 64 x 128). |2(J We actually used data with system 
sizes, TV = 599 x 600 ~ 6399 x 6400 in this paper, and 
checked the finite-size effect and the range of observing 
time in which the system behaves as an infinite system 
at each temperature. Thus, our data can be considered 
as closed to the thermodynamic value. This is a reason 
why we could detect the KT transition accurately. 

In this paper, we have supposed that the ferromag- 
netic correlation along the chain direction remains finite 
or at least it is critical in the IC phase. Thus, the T c \ we 
have obtained is the temperature at which the ferromag- 
netic correlation becomes critical. Above this tempera- 
ture, each ferromagnetic chain exhibits paramagnetism, 
even though the correlation length is very large near the 
T c i. The present results only exclude the finite IC phase 
of this type. In this study, it is clarified that the domain 
walls penetrate into the system along the chain direction 
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at the lower transition temperature, T C 2- The spin con- 
figuration along the chain direction changes drastically 
from the ferromagnetic ordered state to the paramag- 
netic disordered state at the same temperature within 
limits of error. We can neglect the IC phase where the 
ferromagnetic correlation has been expected to be crit- 
ical. Therefore, the domain walls mostly run straight 
along the chain direction. We consider this is why the 
naive free-fermion approximation of Villain and Bak 
gives the best quantitative agreements with our estimate 
of T C 2, and is the best approximation. The analyses for 
other regions of k and determination of the critical expo- 
nents will be a task in the future. The NER of the Binder 
parameter, the specific heat and the spin correlation with 
high accuracy is necessary. Q 

As for the dielectrics, Pb(Zri_ x Ti x )03, the phase di- 
agram Q in the low concentration of Ti is very similar 
to that of the ANNNI model if we interpret the paraelec- 
tric, the ferroelectric and the antiferroelectric phases as 
the paramagnetic, the ferromagnetic and the antiphase 
phases. There are two controversial explanations for the 
existence of the incommensurate phase observed in the 
experiment. Ricote et al. Q| concluded that the appear- 
ance of the IC phase is due to the surface effect by two 
experiments using the powder neutron diffraction and a 
transmission electron microscope. On the other hand, 
Viehland et al. Q considered it a bulk effect by directly 
observing the high resolution image of a transmission 
electron microscope. In addition Watanabe et al. || also 
described that it is not because of the surface effect by 
examining the stability of the IC phase in the bulk. If 
this compound can be explained by the ANNNI model, 
the existence of the IC phase is only possible in three di- 
mensions. Therefore, we predict that the appearance of 
the IC phase in Pb(Zri_ x Ti x )03 is the effect of a bulk. 

In this work, we have presented that the NER anal- 
ysis is very effective for systems with the KT transition 
and/or with slow-dynamics which has been difficult by 
numerical methods. Furthermore, a phase transition is 
studied using the finite-time scaling if we know a proper 
quantity that can probe characteristic time, r, or a char- 
acteristic length, £, even though it is not a relevant order 
parameter. It is considered that images of the phase tran- 
sition which has been believed as standard may change 
in some systems. 
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FIG. 1. The finite-size effect of the layer magnetization for 
k — 0.8 at T = 1.4. Length along an axial direction is fixed to 
800, while that along a chain direction is varied from 799 to 
25599. The converging value is roughly considered as £ X /L X , 
from which we derived the correlation length along the ferro- 
magnetic chain direction as £ x ~ 500. 
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FIG. 3. (a)The finite-time scaling of the antiphase magne- 
tization for k = 0.6, where A = 0.015. We use the data at five 
temperatures ranging T = 0.98 ~ 1.06. (b)The least-square 
fitting of the relaxation time, T( 2 > (e), by Eq. where we obtain 
T c2 = 0.89 ± 0.02, zv = 4.07. 




t[MCS] 

FIG. 2. The NER of the antiphase magnetization, m^ 2 )(t), 
for k = 0.6. The system size is N = 6399 x 6400. 
The antiphase magnetization clearly decays exponentially at 
T > 0.98. 




10 
t[MCS] 

FIG. 4. The NER of the layer magnetization, mi(t), for 
k = 0.6. The size of system is N = 6399 x 6400. The layer 
magnetization clearly decays exponentially at T > 0.98. 
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FIG. 5. (a)The finite-time scaling of the layer magnetiza- 
tion for k — 0.6, where A = 0.015. The data of five temper- 
atures ranging T = 0.98 ~ 1.06 are plotted together. (b)The 
least-square fitting of the relaxation time, Ti(e), by Eq. (where 
we obtain T cl = 0.89 ± 0.02, B = 3.44. 
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FIG. 6. The NER of layer magnetization, m;(t), at 
T = 0.92 and k = 0.6. Three curves of N = 1999 x 2000, 
N = 1599 x 1600 and N = 1199 x 1200 fall on the same line 
until 25000 MCS. Then, the curve of N = 1199 x 1200 deviates 
down from the other two curves, which remain consistent un- 
til 100000 MCS. We discard the data of N = 1199 x 1200 after 
25000 MCS, though we employ the data of N = 1599 x 1600 
until 100000 MCS. 
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FIG. 7. The NER of layer magnetization, mi(t), for 
(a)K = 0.6 and (b)« = 0.8 started from the paramagnetic 
state. 




FIG. 8. The local exponent of the layer magnetization, 
mi(t), for (a)fc = 0.6 and (b)n = 0.8 started from the para- 
magnetic state. In both systems, local exponents converge to 
A = (2 — rf)/z = 0.49 in the limit of t — > oo at the transition 
temperature. 
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FIG. 9. The finite-time scaling of the layer magnetization, 
mi(t), started from the paramagnetic state. The data of four 
different temperatures ranging T = 0.92 ~ 1.00 are plotted 
together for (a)« = 0.6, while those of different three temper- 
atures ranging T — 1.35 ~ 1.50 are plotted for (b)« = 0.8 
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FIG. 10. The least-square fitting of the relaxation time, 
Ti(e), by Eq. ((a)« = 0.6 and (b)/t = 0.8. The transition tem- 
perature is estimated to give the best fitting as (a)T c i = 0.89 
and (b)T c i = 1.30. 
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Present results / References (year) 


k = 0.6 

Td 


T c2 


k = 0.8 


r c2 


NER from T = (scaling) 


0.89(2) 


0.89(2) 


1.31(2) 


1.32(2) 


NER from T — oo (local exponent) 
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1.325(25) 




NER from T = oo (scaling) 
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(1987) 
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.16(4) 



~ 0.91 
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~ 1.10 
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~ 1.05 
0.91(1) 
~ 0.91 



1.45 
1.70 

1.60 
1.65 
1.95 
1.60 



1.30 
1.50 
1.50 

1.35 

1.35 



TABLE I. Comparison the present estimate of T c i and T C 2 for k = 0.6 and k = 0.8 with the previous ones. 
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(from T — 
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1.04 
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1.48 


1.000 


1.000 




1.02 


2.950 


3.025 




1.46 


1.713 


1.715 




1.00 


5.889 


6.145 


1.000 


1.44 


3.033 


3.140 




0.98 


13.301 


16.657 




1.42 


6.264 


6.889 




0.97 






3.202 


1.40 


14.822 


20.190 


17.503 


0.95 






11.103 


1.35 






752.314 


0.92 






500.115 











TABLE II. The relation between the temperatures, T, and the relaxation time, T/2\{e), of the antiphase magnetization, 
m(2) (i), and the relaxation time, Ti(e), of the layer magnetization, mi(t), for k = 0.6 and 0.8. The relaxation time is scaled so 
that a value at the highest temperature becomes unity. 





T 


L X X Ly 


Crossover MCS 


k = 0.6 


0.88 


999 x 1000 


~ 8000 




0.90 


999 x 1000 


~ 10000 




0.92 


1599 x 1600 


~ 100000 




0.95 


999 x 1000 


~ 10000 




0.97 


799 x 800 


~ 8000 




1.00 


999 x 1000 


~ 100000 


k = 0.8 


1.30 


799 x 800 


~ 5000 




1.33 


999 x 1000 


~ 10000 




1.35 


1199 x 1200 


~ 100000 




1.40 


599 x 600 


~ 8000 




1.50 


599 x 600 


~ 100000 


TABLE III. 


The crossover MCS at each temperature for k 


= 0.6 and 0.8. This is 


a MC step until which the NER function 


falls on the same line as that of a larger system. 
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